Fractional Order Hybrid Systems and Their Stability 



S. Hassan HosseinNia, Ines Tejado and Bias M. Vinagre 
Department of Electrical, Electronic and Automation Engineering 
Industrial Engineering School, University of Extremadura, 06006 Badajoz, Spain 
e-mail: {hoseinnia;itejbal;bvinagre} @ unex.es 



Abstract — This paper deals with hybrid systems (HS) with 
fractional order dynamics and their stability. The stability of 
two particular types of fractional order hybrid systems (FOHS), 
i.e., switching and reset control systems, is studied. Common 
Lyapunov method, as well as its frequency domain equivalence, 
are generalized for the former systems and, for the latter. Re- 
condition is used -frequency domain equivalence of Lyapunov- 
like method for reset control systems. The applicability and 
efficiency of the proposed methods are shown by some illustrative 
examples. 

Index Terms — Fractional order hybrid system. Switching 
system. Reset control system. Common Lyapunov method, 
Lyapunov-like method, H^g -condition. 

I. Introduction 

Hybrid systems (HS) are heterogeneous dynamic systems 
whose behaviour is determined by interacting continuous- 
variable and discrete-event dynamics, and they arise from 
the use of finite- state logic to govern continuous physical 
processes or from topological and networks constraints in- 
teracting with continuous control |T|, f2l, O, (H. Typically, 
their stability is analyzed by Lyapunov's theory (see e.g. 151, 
O, Q, O. However, recently a frequency domain method 
equivalent to the common Lyapunov was proposed in |9| to 
analyze the stability of a particular class of HS. 

This paper deals with HS with fractional order dynamics 
and summarizes our stability results proposed in ifTOl , ifTTIl , 
|[T2l for two types of HS: fractional order switching and 
reset control systems. In |[T3l , f\M the applied the proposed 
stability analysis to design a robust fractional-order controller. 
The motivation for studying switching systems came partly 
from the fact that such systems have numerous applications 
in control of mechanical systems, process control, automotive 
industry, power systems, traffic control, and so on, as well 
as there exists a large group of nonlinear systems which 
can be stabilized by switching control schemes, but not by 
any continuous static state feedback control law |[T5]| , |[T6l . 
Likewise, reset control systems arise from overcoming the 
limitations in linear systems. On the other hand, fractional 
dynamics can be found in control systems due to both the 
system itself and the used control strategy (see e.g. ifTTIl . ifTSl ). 
As a result, studying fractional dynamics of HS may be also 
an interesting topic. 

The rest of the paper is organized as follows. Section |ll| 
gives some basic theorems and definitions concerning stability 



of both fractional and integer order systems. In Section III 



stability conditions for fractional order switching systems are 



reset control systems. Section |V| gives some examples to 
show the applicability and goodness of the developed stability 
theory. Finally, Section |Vl| draws the concluding remarks. 

II. Preliminaries 

This section recalls some basic definitions and theorems 
concerning stability of fractional and integer order systems 
which will be useful to present the results of the following 
sections. 

A. Definitions 

Switching systems are hybrid dynamical systems consisting 
of a family of continuous-time subsystems and a rule that 
orchestrates the switching among them IS, |[T9l . On the other 
hand, reset control systems are a class of HS 1201 which 
include a linear controller which resets some of their states 
to zero when their input is zero or certain non-zero values. A 
fractional order linear time invariant (FO-LTI) system can be 
given by: 

D''x=Ax,xeMP (1) 

where a is the fractional order and the operator denotes 
Riemann-Liouville definition given by (see e.g. 1 18J ): 

1 d fit) 



aD^m = 



■ / ^ 

■ Ja (t- 



-dT 



(2) 



r{l - a) dx Ja (t-T) 

where e and t are the lower and upper bounds of the operation, 
a G M is the order and [.] means the integer part. 

B. Stability of fractional order systems 

Theorem 1 (|21 1). A fractional order system ([7]) with order 
a, < a < 1, is t~^ asymptotically stable if and only if there 
exists a positive definite matrix P G such that 



(3) 



where 



-{-AY 



Theorem 2 (| 21|). A fractional system ([7]) with order a, 1 < 
a < 2, is t~^ asymptotically stable if and only if there exists 
a matrix P = P^ >0, P eR'^ 

(A^P + M)sin0 
{-A^P^PA) COS0 



such that 



(A^P-PA) COS0 
{A^P^PA) sin0 



<0, 



(4) 



established. Section IV addresses stability of fractional order 



where = 

Theorem 3 ( 1221 ). A system given by ^ is quadratically 
stable if and only if there exists a matrix P = P^ > 0, P G M"^", 



such that 



A\P^PAi<Oyi=l,...,L. 



Theorem 4 (f9|). Consider c\{s) and C2{s), two stable poly- 
nomials of order n corresponding to the subsystems x = Aix 
and X = A2X, respectively, then the following statements are 
equivalent: 

^) ctI^ ^^nc^fy positive real (SPR). 

2) |arg(ci(7a)))-arg(c2(7a)))| < f , V (O. 

3) Ai and A2 are quadratically stable, which means that 
3P = P^ >0e such that A^P^PAi < , A^P + 
PA2 < 0. 



C. Stability of integer order switching systems 

Consider a switching system as follows: 

x = Ax,A e co{Ai,...,Al} , (5) 

where co denotes the convex combination and Ai, i= 1, ...,L, 
is the switching subsystem, which can be alternatively written 
as El: 



L L 

-Ax,A = Y, ^Ai^Xi > 0, ^ A, = 1 . 

i=l i=l 



(6) 



Theorem 5 (Lyapunov-like theorem |4|). Consider a closed- 
loop reset system given by p2] ). If there exists a Lyapunov- 
function candidate V{x) such that 



V{x) < 0, x{t) ^ . 



(7) 



where Ap G R"p><"^ Bp G 



plxn, 



^ Ac e 



Be eW''''^ and Q G R^'''''. Therefore, the closed-loop reset 
control system can be then described by the following FDI: 

D^x{t) =Acix{t) ^Bcir, x(t) ^ 

x{t^) =ARx{t), x{t) G ^ (12) 
y{t) = Ccix(t) 



where t^ denotes t and k is sample time, x = 



Ad 



BpCc 



"4. 





Ar^ 



Ccl ■ 



IS 





Ac BcCr 

BrCp A, _ 
\Cp 0] and Bd = [O By^ - The reset surface 
defined by: 

^ = {xeW \ Ccix = {), {I-Ar)x^Q}, (13) 
where n = ny^nc^np. 

III. Stability of Fractional Order Switching 
Systems 

Our objective hereafter is to establish stability conditions for 
fractional order switching systems. In this section, we firstly 
present the stability of such systems by common Lyapunov 
functions, which have been previously generalized to frac- 
tional order switching systems, and further its equivalence in 
frequency domain. This developed theory can be found in |[Tq1 , 

mi. 



Ay(x) = y(x(^+))-y(jc(0) <o, x{t) g^, 

then there exists a left-continuous function x{t) satisfying ( [72] ) 
for all ^ > 0, and the equilibrium point x^ is globally uniformly 
asymptotically stable. 

D. Dynamics of reset control systems 

The dynamics of a reset controller can be described by a 
FDI equation as: 

D'^Xr(t)=ArXr(t)^Bre{t), e(t)^0, 

Xr{t+)=ARMt). e{t)=Q, (9) 

Ur{t) =CrXr(t), 

where < a < 1 is the order of differentiation, Xr(t) G M"^ 
is the reset controller state and Ur{t) G M is its output. The 
matrix Ar^ ^ j^^rx^r identifies that subset of states Xr that are 

reset (the last ^ states) and has the form Ar^ — ^ 



(8) A. Common Lyapunov theory 

Consider a fractional order switching system of the form of 
([5]) as 

D^jc = Ax,A G CO {Ai , A^} . (14) 

Theorem 6. A fractional system described by ( fM) ) with order 
a, 1 < a < 2, is stable if and only if there exists a matrix 
P = pT >0, Pe M"><^, such that 



0; 



with n^ = nr — n^. The linear controller C{s) and plant P{s) 
have, respectively, state space representations as follows: 



and 



D'^Xc{t) =AcXc(t)^BcUr(t), 
Uc{t) =CcXc{t), 

D'^Xp{t) =ApXp{t)^BpUc{t), 
y{t)=CpXp{t), 



(10) 



(11) 



{aJp + PAi) sin {aJp - PAi) cos 0' 
-AJp + PAi) cos (aJp + PAi) sin 



<0,V/=1,..,L, (15) 



where (j) = 

Theorem 7. A fractional system given by (14^ with order a, 
< a < 1, is stable if and only if there exists a matrix P = 
P^ >0, Pe W"", such that 



P^PM<0, V/: 



1,...,L. 



(16) 



B. Frequency domain approach 

Next, frequency domain stability conditions will be given 
for fractional order switching systems based on results in (91. 
Consider a stable pseudo-polynomial of order na of system 
( [T4| ) as 

d{s) = + + • • • + Ji^^ + Jo, (17) 



and a polynomial of order n of system x = Ax as 



c{s) 



-Cn-lS 



(n-1). 



-ClS + CQ. 



(18) 



In the following, the necessary and sufficient condition for 
the stability for fractional order switching systems will be 
given. 

Theorem 8. Consider di{s) and d2{s), two stable pseudo- 
polynomials of order n corresponding to the subsystems D^x = 
A\x and D^x = A2X with order a, 1 < a < 2, respectively, then 
the following statements are equivalent: 
1) |arg(det((A^-(O^/)-27(oAisin0)) - 

arg(det((A^-(O^/)-2j(oA2sin0))| < f ,V(0, 
being I the identity matrix with proper dimensions. 
A] and A2 are stable, which means that 3P = P^ >0e 



2) 



^1 

such that 



[AjP^PAi) sin0 
{-Aj P^PAi) cos (j) 



{AjP-PAi) cos (I)' 
{AfP^PAi) sin0 



<0,V/=1,2. 



Theorem 9. Consider two stable fractional order subsystems 
D^x = Aix and D^x = A^x with order a, < a < 1, then the 
following statements are equivalent: 

1) |arg(det(M-jW))-arg(det(j2/2-jW))| < f , V to. 

2) A\ and A2 are stable, which means that 3P = P^ > G 

such that 

j^^P + Pj^<0,V/=1,2. 

Although the theory developed in the frequency domain 
does not necessarily prove the SPRness, a relation equivalent 
to the stability was obtained. 

IV. Stability of Fractional Order Reset Control 
Systems 

In this section, fractional reset control systems will be 
handled by fractional order differential inclusion (FDI) equa- 
tions. Stability of this kind of systems will be analyzed using 
Lyapunov-like method presented previously. This developed 
theory can be found in |12|. 



Definition 1. Reset control system {12) is said to satisfy the 
-condition if there exists a j3 G M"^ and a positive-definite 
matrix P^ G ]R"^><"^ such that 



H(^{s)=[l5Cp On^ Pm\{sI- 



' 
1m 



(19) 



where sz/ = (-A^/) 2-« ^. 

According to 1241 . 1251 , 1261 , an integer-order reset control 



system of the form of (12) -with a = 1- is asymptotically 
stable if and only if it satisfies the H^g -condition. The same 
idea can be used to prove the stability of fractional-order reset 
systems. 

Now, consider y{z(t)) = z(tY 3^z(t), ^ e R^""^ as a 
Lyapunov candidate for the unforced reset system ([12]) (r = 0) 
where x= [0, • • • ,0, l]z(0, z{t) G M^^^^, z = Afz{t), and 



" 



aV«- 





(see 1211 more details for this 



ai/« 

transformation). Then, in accordance with f2T\, the necessary 
and sufficient condition to satisfy V{z{t)) <0 when i < a < 1 

is: (Aa^ P^p(Aa^ < 0, x{t) ^ J{ , where P(c ^) G 
R"^" > 0. Likewise, based on results stated in Theorem [H 
the necessary and sufficient condition for < a < 1 is 

j^^P + Pj^<0, x(t) ^Ji. 

Transforming the second equation of reset system ([12]), we 
have 



iN-n 

Ar 



zit), 



(20) 



where If^-n is identity matrix with dimension of N — n. Thus, 

Ay(z(f)) <o if 









_ Al 




Ar 



z{t) < 0. (21) 



Then, ([21| is satisfied if V{x(t"^)) -V{x(t)), 

x^{t){AlPAR-P< 0)x{t) < 0, x{t) G 

Therefore, Theorem [5] can be reshaped in the following re- 
mark. 

Remark 1. Choosing V{z) = z{tY S^z{t), ^ G as a 

Lyapunov candidate, and applying Theorem^ fractional- order 
reset system ( [72] ) is asymptotically stable if and only if: 

(22) 
(23) 



aIpAr-P<0, x{t) G. 



Remark [T] can be also applicable transforming x to x(t) = 
x(t)—Xe = x{t) ^ A~^Bcir in (12) in a special case of the 



reset control system with constant input r. Define ^ — 
{x G : Ccix(t) = r}, and let <I> be a matrix whose columns 
span Since ^ C ( |23| ) is implied by 

c^{aIpAr-P <0)c^<0. 



(24) 



A straightforward computation shows that inequality ( [24] ) 
holds for some positive-definite symmetric matrix P if there 
exists a j3 G and a positive-definite P^ G R"^><"^ such 
that 



[0 0^ I^]P=[l5Cp On^ P^] 



(25) 



To analyze stability, it suffices to find a positive-definite 



symmetric matrix P such that ( [22] ) and ( [25] ) hold. Taking into 
account Kalman-Yakubovich-Popov (KYP) lemma 1271 . such 



P exists if (^) in ( 19 ) is strictly positive real (SPR) for some 
j3. In addition, in accordance with |28|, it is obvious that the 
H^(^) is SPR if 

|arg(if^(7C0))|< J,Vco. (26) 



Therefore, these results can be stated in the following theorem. 

Theorem 10. The closed-loop fractional-order reset control 
system ( [72] ) is asymptotically stable if and only if it satisfies 
the Hj^ -condition ( [79| ) or its phase equivalence ( [26| . 

V. Examples 

This section gives some examples in order to show the ap- 
plicability and effectiveness of the stability theories developed 
for FOHS. To this respect, phase portraits and time responses 
of the systems will be shown. 



Example 1. Consider the switching system ( |74| ) with L 

■1 1 s 11 ■ . T-O.l 0.1 ^ 

with the jollowmg parameters: A\ - 



-2.0 -0.1 



A2 



-0.01 
-0.1 



2.0 
-0.01 



and order a, < a < 1. 



Applying Theorem |9] the phase difference condition should 
be satisfied for all a, < a < 1, to guarantee the stability -this 
condition is depicted in Fig. [T] for < a < 1 with increments 
of 0.1. As can be seen, the fractional order system is stable 
for a G (0,0.6]. The phase differences when a G [0.7,1] are 
bigger than 7t/2 which indicates unknown stability status, i.e., 
the system may be stable or unstable. For better understanding 
of this initial notice on the system stability, its phase portrait 
is shown in Fig. [3] for three values of a (a = 0.6, a = 0.8 
and a = 0.9). The green trajectory is an example to show the 
stability or instability of the switching system. The following 
conclusions can be stated from these results: 

• When a = 0.6, it can be observed that the system is stable 
for arbitrary switching. This can be also confirmed by 

r 1 0.2" 

0.2 1 



satisfies the stability 



-1.47 
-0.89 

-1.45 
0.59 



-0.65 
-0.54 

0.59 
-0.47 



<0, 



<0. 



the fact that a matrix P 
conditions as follows: 

(-(-Al)^)^P + P(-(-Al)r4): 

(-(-A2)^)^P + P(-(-A2)^) 

The switching region is shown in Fig. [2] in which Ci 
refers to the zone which only subsystem 1 is active, 
whereas C2 is the zone which corresponds to subsystem 
2. Z) is a common region with a random layer where both 
system can be active. The red lines indicate the switching 
from subsystem 1 to subsystem 2, whereas the blue lines 
show the switching in contrary. 

When a = 0.8, its phase portrait shows almost the same 
behaviour as with order a = 0.6. Although one cannot 
find a trajectory which leads to unstable switching system, 
the stability of the system under arbitrary switching is on 
doubt. 

Finally, in the case of a = 0.9, the system will be unstable 
if it switches like the green trajectory shown in Fig.jSjc). 



Example 2. Now, 
{Mi with L = 2 

-0.2 -1.0" 

0.01 -0.1 



A2 



consider the switching system given by 

with the following parameters: Ai = 

-0.3 0.01] ^ ^ , 
^ Q Q ^ and order a, 1 < a < 2. 




Fig. 1: Phase differences of characteristic polynomials of 
system in Example [T] for different values of its order a, 
0< a < 1 





\ r 















Fig. 2: Switching region for random switching of system in 
Example [T] 



It is easy to find that the subsystem 1 is stable for when a G 
(1,1.67), whereas the subsystem 2 is stable for all values of 
a G (1,2). Therefore, applying Theorem [s] when a G (1, 1.67), 
the following condition 



arg ( det 



arg det 



0.03-0)^ + 70.40) sin 
-0.003 -70.020) sin 

0.08 -0)2 + 70.60) sin 
0.4 + 72o)sin0 



0.3 + 72o)sin0 
-0)2 + 7O.2o)sin0) 

-0.004 -7O.O20) sin 
-0)2 + 7O.2o)sin0) 



<-,Vo) 
(27) 

should be satisfied, V a, 1 < a < 1.67. The phase difference 
( [27] ) is depicted in Fig. |4ja). In order to make the results 
clearer, the maximum values of ([27]) are also plotted in 
Fig. |4jb) versus the order of the system. It can be seen that the 
system is stable if a G (1,1.65). The stability of the system 
when a G [1.65,1.67) is unknown. 

Example 3. Let us consider the same feedback system as in 
h29] with the system, the base controller and reset controller 



{a) 






Fig. 3: Phase portrait of system in Example [T] when: (a) a = 
0.6 (b) a = 0.8 (c) a = 0.9. The blue trajectories refer to 
subsystem 1 , whereas the red ones correspond to subsystem 2 



C{s) 



1 and, R{s) 



transfer functions P{s) = ^2^q 
^oT^y respectively. Stability analysis of CI (b = O^a = I), 
FORE (b y^O, a = 1) and FCI (b = 0,a = 0.5) will analyzed 
in this example. 

The system stability will be analyzed for the first order reset 
element (FORE) controller, the Clegg integrator (CI) and the 
fractional Clegg integrator (FCI). For FORE controller, the 
integer-order closed-loop system can be given by: 



=ArX -- 



1 
-0.2 1 x{t) 
-1 -1 -b 

1 Ol 

1 x{t) 



y = C,ix=[f 1 0]x{t) 



1 1<a<1.65 


120 




^ 1 a=[1 .65, 1 .66] 
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Frequency (rad/sj 



1 1.1 1.2 1.3 1.4 1.5 1.6 



Fig. 4: Stability of the system in Example[2]for different values 
of its order a, 1 < a < 2: (a) Phase difference of condition 



( 27 ) (b) Maximum value of ( 27 ) versus a 



where x(t) = [xp^{t)^Xp2{t),Xr{t)]^ . And, the closed-loop sys- 
tem using FCI can be stated as 



,0.5 , 
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1 





0" 
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-0.2 


1 




-1 





-1 


0_ 




) = 


' h 

0l,4 


04,r 






y = C,i^{t)=[\ 





1 








where ^{t) = [^p,{t)r ■ ■ . ^p,{t)Mt)Z^ ^Piit) = ^Pii^)^ 



^p^{t) =Xp^{t). According to condition (19), Hp correspond- 
ing to FORE and FCI are simply given by (Tor both case FORE 
and FCI n^ = 1 and then = 1): 



[P I]{sl-A, 



52 + 0.25 + 0.8j8 



' + + 0.2)52 + (1+0.2Z?)5+1' 



and 



Hl^\s)=[P /3 l](./-(-(-A,,)i))" 



(28) 



(29) 



Therefore, using Theorem 10 
trolled by FORE and FCI 

jjFORE 
H/3 



the closed-loop systems con- 
are asymptotically stable if 
28l), 



[s) and H^^^(^) are SPR. Substituting /? = 1 in 
the FORE reset system is asymptotically stable for all 0.42^ 
J3 < 1.46. With respect to CI (similarly to FORE with b = 0), 
stability cannot be guaranteed with this theorem. And applying 
FCI, it can be easily stated that the system is asymptotically 
stable for j3 < 0.62. In addition, the phase equivalences 
corresponding to (28 ) and (29) are shown in Fig. [5] for j3 = 0.5 
and b = I. It can be seen that both phases verifies condition 
(|26l), which has concordance with the theoretical results. 




10"^ 10"^ 10° 10^ 10^ 

Frequency (rad/s) 



Fig. 5: Phase equivalence of H^g in Example pi (a) Applying 
FCI (b) Applying FORE 



VI. Conclusions 

This paper has addressed stability for two classes of hybrid 
systems (HS), i.e., switching and reset control systems, with 
fractional order dynamics by means of a Lyapunov like 
method. More precisely, common Lyapunov method as well 
as its frequency domain equivalence were generalized for the 
fractional order switching systems. Likewise, H^g -condition 
was used -frequency domain equivalence of Lyapunov like 
method for reset control systems- to prove stability for 
fractional order reset control systems. The applicability and 
efficiency of the proposed methods were shown by some 
illustrative examples. 
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